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A CENTRAL LIMIT THEOREM FOR FLUCTUATIONS IN ONE 
DIMENSIONAL STOCHASTIC HOMOGENIZATION 

YU GU 


Abstract. In this paper, we analyze the random fluctuations in a one di¬ 
mensional stochastic homogenization problem and prove a central limit result, 
i.e., the first order fluctuations can be described by a Gaussian process that 
solves an SPDE with additive spatial white noise. Using a probabilistic ap¬ 
proach, we obtain a precise error decomposition up to the first order, which 
helps to decompose the limiting Gaussian process, with one of the components 
corresponding to the corrector obtained by a formal two scale expansion. 
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1. Introduction 
The equation we are interested in is 

1 X 

(1.1) d t u e = -d x a(-,uj)d x u e , t> 0,ieK, 

2 e 

with an initial condition u e (0,x) = f(x) G C£°(R). Here a is a smooth stationary 
random field defined on some probability space (H, J 7 , P), and satisfies 

A < d(x, oj) < 1 

for some A £ (0,1) and i£l,we!l. The standard homogenization result shows 
that u e Whom with 

(1.2) C^Uhom = 

and the homogenization constant a is the harmonic mean of a: 

a -1 = E{a -1 }. 

The goal is to analyze the first order fluctuations, i.e., if the homogenization is 
viewed as a law of large numbers type result, we are interested in a central limit 
theorem (CLT) here. The same question has been addressed for the boundary value 
problem 

(1.3) -= f, with a; € (0,1), 

and u £ (0) = 0,u e (l) = 1, under different assumptions on the correlation properties 
of a [4, 2, 9]. It was shown in [4, Theorem 3.1] that if a satisfies certain mixing 
assumption, a CLT holds: 

u e (x) - u hom (x) ^ r 1 F ^ y)dBy 
V e Jo 

in C([0,1]), where F(x, y ) is deterministic and B y is a standard Brownian motion. 
The analysis used the explicit solution to (1.3) and reduced the problem to the 
weak convergence of oscillatory random integrals. We “revisit” the problem on the 
whole space for the following reasons: 


l 
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(i) It seems the approaches in [4, 2, 9] fails to work for (1.1) since we do not have 
an explicit solution when the problem is posed on the whole space, so a different 
method needs to be developed. 

(ii) It is not clear how the boundary condition affects the asymptotic behavior 
of the rescaled fluctuations in (1.3), and the Wiener integral f Q F[x,y)dB v is not 
naturally related to the corrector obtained through the formal two scale expansion 
that is used extensively in homogenization. In this paper we are looking for an 
error decomposition that is in parallel to the formal expansion and also indicates 
clearly how each term contributes to the limiting Gaussian distribution. 

(iii) It was shown recently when d > 3 that a pointwise formal two scale expansion 
holds [11, Theorem 2.3]: for fixed ( t,x ), 

~ X 

(1.4) u e {x) = Uhom{x) + eVuhom(i) ■ </>(-) + o(e), 

£ 

where (f) is the mean-zero stationary corrector and o{e)/e —> 0 in L 1 (Cl). It can be 
seen from (1.4) that the local (pointwise) fluctuation is not necessarily Gaussian 
since </> is not Gaussian. On the other hand, the large scale central limit results are 
proved for the global fluctuations of the solutions and correctors [12, 17] 

(1.5) -^72 [ (u e {x) - E{u e {x)})g(x)dx =>• lV(0,crf), 

£ ' JR d 

(1.6) [ e Vu ho m(a;) ■ <t>(-)g(x)dx => iV(0,of)- 

It turns out that o\ ^ <72, so the corrector which represents the local fluctuation 
does not capture the global fluctuation! Mathematically it is not surprising since 
the o(e) in (1.4) could contribute on the level of e d / 2 when d > 3. From a practical 
point of view, it is important that we can extract the right term from o(e ) such that 
together with the corrector they represent the global fluctuations of the solutions. 
To understand the mechanism better, we start from the simpler setting d = 1, where 
the local and global fluctuations are known to be described by a single Gaussian 
field on the level of yfe. We expect the error decomposition and its relation to the 
corrector in low dimension to shed light on the situation in high dimensions. 

Quantitative aspects of stochastic homogenization of divergence form operator 
has witnessed a lot of progress recently, from both analytic and probabilistic points 
of view [7, 5, 6, 8, 16]. Our approach falls into the more probabilistic side: we 
use the probabilistic representation of the solutions to (1.1) and quantify the weak 
convergence of an underlying diffusion in random environment. The main ingre¬ 
dients in our analysis consist of the Kipnis-Varadhan’s method [14, 15] applied to 
reversible diffusion in random environment and the quantitative martingale CLT 
developed in [16, 10] to extract the first order error in the martingale convergence. 
We also rely heavily on the work [13], where the authors analyzed the asymptotics 
of 

(1.7) dtu £ — —D x d(—, uj'jdx'Uz -f- ^c( , w)i£ e , 

i.e., (1.1) with a large highly oscillatory random potential. It turns out that a part 
of the error in our analysis of u e — Whom solves (1.7) with an additive rather than 
multiplicative potential. By following a similar argument, we obtain a limiting 
SPDE (for this part of the error) with additive white noise (which is a Gaussian 
process), while the limit of (1.7) is an equation with multiplicative white noise. 

The Kipnis-Varadhan’s method decomposes the underlying diffusion process as 
a small remainder plus a martingale which converges to the limit. One of our 
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main contributions is to combine the errors coming from the remainder and the 
martingale convergence, which are three Gaussian processes, and write the sum as 
the solution to one SPDE with additive noise. On one hand, this justifies rigorously 
in the simpler setting d = 1 the heuristics presented in [12]. More precisely, if we 
assume that in (1.1), the fluctuations of the coefficient a around its homogenization 
limit a can be described by some mean zero, strongly mixing processes, denoted by 
V in (2.1) below, and on the large scale, V can be replaced by some spatial white 
noise VV (after proper rescaling), then (1.1) may be rewritten as 

(1-8) d t u e ft ^ d x {a + \feVV)d x u e , 

so the error v s = £ _1 / 2 (u e — «hom) should solve 

(1-9) d t v e ft ^ad x v e + ^d x (d x u hom W). 

Indeed, we obtain the limiting fluctuation described by an equation of the form 
(1.9) in Theorem 2.2 below. For more detailed discussions of the heuristics in high 
dimensions, we refer to [12, Section 1]. On the other hand, the three Gaussian pro¬ 
cesses obtained in the limit provide a natural decomposition of the limiting SPDE 
(on the level of equations), which is in parallel to the martingale decomposition 
of the underlying diffusion in random environment (on the level of stochastic pro¬ 
cesses), and this helps us to see the role played by the corrector more clearly. It is 
not clear whether (1.9) can be obtained through some PDE arguments. 

Here are some notations used throughout the paper. The expectation in (f2, J r , P) 
is denoted by E. When averaging with respect to some independent Brownian 
motion B 1 W, we use E^, Ew- The normal distribution with mean fi and variance 
a 2 is denoted by N(n, <r 2 ), and the density function of 7V(0, t) is denoted by q t (x ) = 
(-\/27rf) -1 e - l x l / 2t . We write a < b when a < Cb for some constant C independent 
of t. x , £. 

The rest of the paper is organized as follows. We present the setup and main 
results in Section 2. The error decomposition is discussed in Section 3, and weak 
convergence results are obtained in Section 4. In Section 5 we finish the proof of the 
main result and compare it with high dimensions d > 3. Some technical lemmas 
are left in Section A. 


2. Setup and main results 

We first assume there is a group of measure-preserving, ergodic transformation 
{t x , x € R} associated with the probability space (f2, J-, P), then the coefficient field 
a(x,to) is defined by 

a(x, to) = a(T x to) 

for some a : 0 —> [A, 1]. We further assume it is smooth and of finite range 
dependence: 

(i) a(x, to) has C 2 sample paths whose first and second order derivatives are 
uniformly bounded in (x,oo). 

(ii) For any two sets 4,BC1, if dist(A,H) > 1, then Ta = a(d(x,to) : x £ A) 
is independent of Tn = a(a(x,tu) : x € B). 

Remark 2.1. The finite range of dependence can be replaced by some mixing con¬ 
dition, e.g., the (j )—mixing used in [13]. 
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Besides the coefficient field a(x,u;), another important random field in our anal¬ 
ysis is 


( 2 . 1 ) 


V(x,u) 


1 1 

a(x,ui) a’ 


which may be seen as the fluctuations of the homogenization constant. It is clear 
that V is of finite range dependence, and its covariance function is given by R(x ) = 
E{F(0, ui)V{x, w)} and the power spectrum is 


( 2 . 2 ) 



R(x)e l ^ x dx. 


The following is our main result: 


Theorem 2.2. Let v e = e 1 ' 2 {u e — Uhom) and v solves 

(2.3) d t v = ^ ad^v - ^R{0)^a 2 d x (d x u hom VV), with v(Q,x) = 0, 

where W is spatial white noise, then as e — > 0, v e => v in the following sense: 

(i) as a process in (t,x), the finite dimensional distributions converge, 

(ii) for any test function g £ C£°(R), f R v e (t, x)g(x)dx =£* f R v(t,x)g(x)dx in 
distribution. 


It turns out the Gaussian process v(t, x) is a superposition of three Gaussian 
processes, and one of them takes the form d x u\, om {t, x)W(x), which corresponds to 
the corrector obtained through a formal two scale expansion. Here 

W{x) := f W(y)dy 
Jo 

is a two-sided Brownian motion. 


2.1. Diffusion in random environment. Our starting point to prove Theo¬ 
rem 2.2 is a probabilistic representation. For every fixed u> £ Ll,x € R and e > 0, 
we define the underlying diffusion in random environment by the Ito’s SDE: 

(2.4) dX t = b(X t ,uj)dt + d(X t ,uj)dBt, with Xq = x/e, 

where 

b{x,u>) = —d'(x,u>), d(x,to) = a^(x,u>). 

The driving Brownian motion B t is built on another probability space (S,^4, Ps). 
It is straightforward to check that eX t / s 2 is a Markov process starting from x with 
the generator L u = i d x a(x/e,u))d x , so the solution to (1.1) can be written as 

u e (t,x) = E B {f(£X t/e 2 )}, 

where Es denotes the expectation in (E,_4 ,P_b). 

It can be shown that eX t / e 2 converges in distribution to x + dWt for some 
Brownian motion Wt starting from the origin and d = \fd, so 

u e (t,x) = E B {f(eX t/e 2 )} E w {f(x + dW t )} = u honi (t,x) 

in probability. It is clear that to further get the first order fluctuations v e = 
£ _1 ' 2 (ii e — Whom), we need to quantify the weak convergence of eX t / e 2 => x + dWt 
up to the first order. 

We define an environmental process by u>t = t_x t w, and it satisfies the following 
properties [15, Proposition 9.8]: 

Proposition 2.3. (w s ) s >o is a Markov process that is reversible and ergodic with 
respect to the measure P. 
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In Sections 3 and 4, we will only show that for fixed (t,x), v E (t,x) => v(t,x) in 
distribution. To simplify the presentation, we will shift the random environment w 
by x/e without changing the distribution of u E (t,x ) (for fixed (t,x)), thus in the 
following we will assume 

u e (t, x) = E B {f{x + £X t/s 2 )} 

with X t solving (2.4) but starting from the origin: 

(2.5) dX t = b(X t ,uj)dt + a(X t ,uj)dBt, with X 0 = 0. 

The convergence of finite dimensional distributions and the global weak convergence 
are discussed in Section 5. 

To simplify the notations, we will omit the dependence on w from now on. 


3. Quenched invariance principle and error decomposition 


To quantify the weak convergence, we first present a proof of eX t / e 2 => dWt , 
where the diffusion in random environment is decomposed as a remainder plus a 
martingale, and the speed of weak convergence hinges on how small the remainder 
is and how “close” in distribution the martingale is to the limiting Brownian motion. 

By the uniform ellipticity condition, we have the following standard heat kernel 
estimates [19] which will be used extensively in our analysis. 

(i) The density function q E (t,x) of eX t / e 2 satisfies 

( 3 . 1 ) qe ^ x )<l. e -o\^/t 

uniformly in £,w for some c > 0. 

(ii) For any t > 0, 

(3.2) P B { sup \eX s/e 2 \ >M}< e ~ cM2/t 

se[Q,t] 

uniformly in e,w for some c > 0. 

The following result is classical. For the sake of convenience, we present the 
proof here. 


Proposition 3.1. For almost every realization u> € Cl, eX t / e 2 => aWt 


Proof of Proposition 3.1. First we write the SDE as the integral equation 


/■t/e 2 


eX^j e 2 e 

•to 

and by solving a corrector equation 
(3.3) 

and applying Ito’s formula, we have 


b{X s )ds + e 


i/e 2 


d(X s )dB s 


1 d_d 

2 dx dx 


l 1 d ~ 


(3.4) £ X t/E 2=RI + Mf, 

with 

R e t = —£(j){X t / E 2 ) + £(j)(X 0 ), 

and 

/•t/e 2 l 

M t =a£ dB s . 

Jo &\X S ) 

Here (f> satisfies 

4>'{x) = d - - = aV(x), 

\a{x) aj 
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and we choose 
(3.5) 


4>{x) = a V(y)dy. 
Jo 


For the remainder, we have 


sup \Rt \ < sup e\(j){£X t / e 2 /e)\ 
te[o,T] te[o,T] 

— SU P e l < ^’( £ ^t/e 2 / £ )l lsup te r 0 T , \eX t/ 2 \<M 

te[o,T] 

+ Sup ^ ^l^(^-^'t/e 2 /£)|l SU p tg j 0 \eX t / e 2\>M 


for any constant M > 0. By ergodicity and the fact that E{V} = 0, we have for 
almost every realization u e f 1, (j>{x)/x — > 0 as |x| — > oo, i.e., the corrector has a 
sublinear growth. This implies 

Sup £\4>(£X t/s 2/£)\l sup \ e x t . 2 \<M 0 

te[o,r] 

as e —> 0. For the second part, we use the sublinear growth to obtain 

£ |(^( e ^Q/e 2 / e )|lsup tg [ 0|T ] \eX t/e2 \>M \^X t / e 2 | l S up tg[0 T] \eX t/s2 \>Mi 
so we only need to show that 


Eb{ sup ^ |£^ft/£ 2 |lsup tg[0 T] 


\eX t/e2 \>m} 0 


uniformly in £ as M —> oo, but this comes from (3.2). To summarize, we have 
shown that for almost every ui, the remainder Rf —> 0 in C(R+). 

For the martingale part, the quadratic variation can be written as 


2 Jl 


r t/e 2 


(M e ) t = a 2 e 


1 


-ds. 


Jo a(r Xs w) 

so by ergodicity of tx s x> (it is independent of e since we have shifted by environ¬ 
ment), we have that for almost every ui, (M e ) t —> at almost surely as £ —> 0, and 
this implies M t e => dWt in C(R+). The proof is complete. 


Now we can decompose the error in homogenization according to the martingale 
decomposition of eX t / e 2 . By (3.4), we write 

u E {t , a:) - u hom {t, x) = E B {f(x + R E t + M t £ )} - E w {f{x + dW t )} := £\ + £ 2 
with 

£i = E B {f(x + Rt + M £ )} - E B {f{x + M £ )}, 

and 

£ 2 = E B {f(x + M e t )} - E w {f(x + dW t )}. 

The main contribution from £\ can be extracted by a Taylor expansion, i.e., we 
have 

\£i - E B {f\x + Mt)R E t }\ < E B {\Rl\ 2 }. 

The main contribution from £ 2 can be extracted by a quantitative martingale CLT 
[10, Proposition 3.2]: 

\£ 2 - ^E B {f"{ x + Mt)((M e ) t - at )}| <E s {|(M £ ) t - dt\i} 

<(E B {\(M E ) t -dt\ 2 })i. 

Now we define 

Vi ,e(t,x) =E B {f'(x + M^)Rt}, v 2 ,e{t,x) = |e B {f"(x + Mf){(M E ) t - at)}, 
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By Lemmas 3.2 and 3.3 below, we have 

(3.6) E{| U £ (t,x) - Uhom(t,x) - VI,e(t,x) - V 2 ,e(t,x)\} < yfl, 

therefore, to analyze the asymptotic distribution of £~ x l 2 [u e — Whom), we only need 
to consider ui )£ + v 2:E . 

Lemma 3.2. EE B {|f?f | 2 } < £\/t 

Lemma 3.3. EE B {|(M e )t — at| 2 } < eti 

Proof of Lemma 3.2. Since = — £(p(X t / e 2 ) + £</>(Xo) and 4>{x) = a f/ V(y)dy, we 
only need to consider £(j)(X t / e 2 ). If we write 

Eb{ \ef)(X t/e 2 )\ 2 } = f \ef)(-)\ 2 q e (t,x)dx, 

Jr £ 

with q e (t,x) the density function of eX t / e 2 (which depends on the random realiza¬ 
tion w), by the heat kernel bound (3.1) we have 

E B {\sf>(X t/e2 )\ 2 } < [ | 

Jr e v t 

for every u>. By Lemma A.l, E{|0(a;)| 2 } < |x|, so we can take E on both sides of 
the above expression to obtain 

EE B {|i?f| 2 }< £ [ |ir|-U- c|x l 2/t dx<eVt. 

Jr v t 

□ 


Proof of Lemma 3.3. First, we write 

rt/e' 


1 


1 


rt/e 2 


ds = a 2 e 2 / V(X s )ds, 
Jo 


(3.7) {. M s ) t -dt = d 2 e 2 . 

Jo \a{X s ) a , 

where we recall V = d _1 — a -1 . Since V has mean zero, by ergodic theorem, we 
have ( M e ) t — at —> 0 as e —> 0, but to quantify how small it is, we need to apply a 
martingale decomposition again, in the same spirit as for b. 

Let if satisfy 


(3.8) 

then 


1 d . d ~ ~ 

-~2ii a di* = V - 


rt/e 2 


V(X s )ds = -e 2 'ijj{X t / e 2 ) +e 2 Jjj(X 0 ) +£ 2 


rt/e 2 


${X a )5{X,)dB s . 


Since 4>{x) = a [/ V(y)dy , we can choose 


(3.9) 


4>(x) = 


4>{y) 

a(y) 


dy. 


By Lemma A.l we have E{| , 0(x)| 2 } < |x] 3 , and we follow the same discussion as in 
the proof of Lemma 3.2. For £ 2 %f(X t / e 2 ) we have 

E B {|e 2 ^(X t/£ 2 )| 2 } = f £ i \^(-)\ 2 q e {t,x)dx, 

Jr £ 

with q e (t , x) the density of £X t / e 2 . Applying the heat kernel bound (3.1) and taking 
E, we conclude that 

EE B {|e 2 i/!(X t / e 2 )| 2 } < £ f \x\ 3 -^=e~ c \ x \/*dx < £t%. 

Jr v t 
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Since Xo = 0, we have £ 2 i/)(X o) = 0. For the martingale term, we have 
/•t/e 2 /-t/e 2 

E s {|e 2 / $(X s )d(X s )dB s | 2 } =e 4 / E B {|^'(X s )d(X s )| 2 }ds 

Jo Jo 

rt/e 2 

<£ 4 / E B {|</>(X s )| 2 }ds. 

Jo 

By a similar discussion, we have 

pt/e 2 nt 

£ 4 / EE s {|^(X s )| 2 }ds =£ 2 / EE s {|^(£X s/e 2/£)| 2 }ds 


<£ 2 / / MJ_ e - C M 2 / 5 dxds<£tJ, 

£ V s 


which completes the proof. 


□ 


4. Weak convergence 

Now we consider Wi l£ ,U 2 ,£ and prove the weak convergence of e~ 1 ' 2 (vi iS + u 2 , e ). 
Let us recall that 


(4.1) v lt£ {t,x) = E B {f\x + Mt)R £ t } 

(4.2) v 2 , e (t, a) = ^{/"(x + M t e )((M E ) t - at)} 

and give a heuristic explanation of what we may expect from the weak convergence 
of 

For fi lE , we can write 

^ = -V~ej>(X t/£2 ) = -^( f ^). 

V £ e 

On one hand, since 4>{x) = a f* V ( y)dy and V is a mean-zero stationary random 
process with finite range of dependence, by a classical functional central limit the¬ 
orem [3, pages 178,179] we have 

y/E<f>( — ) =4> cW(x) 

£ 

weakly in C(R), where 

(4.3) c=R(0)h, 


and W(x) is a two-sided Brownian motion with W(0) = 0, i.e., 



Wi (x) 
W 2 (-ar) 


x > 0, 
x < 0, 


where Wi, W 2 are independent Brownian motions starting from the origin. On the 
other hand, by Proposition 3.1, we have 


eX t / e 2 => a\V t 

in C(R+) for almost every u. Apparently W and W are independent since they live 
in (f2, J 7 , P) and (E,A,Ps) respectively, so we expect 

^ 7 = = -Ve4 >{——) => -cW(dW t ) 

V £ £ 

in distribution. 

For v 2 ,e, we have 


{M e ) t - at 





e 2 


£ 


)ds, 
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and since 

4= [ V{ V -)dy^m?W{x), 

V £ Jo £ 

we can formally write 

y¥)^(0)^ Wl 

V £ e 

where VV is the spatial white noise, and this implies 

/ y( ° s/e2 )ds^ R(0)* / W(aW s )ds = £(0)* / L t (x)W(dx), 

v £ Jo £ Jo Jr 

with L t (x ) the local time of aWt and VC(dx) interpreted as the Wiener integral. 

The above heuristic argument has already been made rigorous in [13, Theorem 
3.1]. More precisely, it was shown that for fixed t > 0, 

(eX t/s 2 ,^= [ V(—^)ds)^(dW u R(0)i [ L t (x)W(dx)), 
v £ Jo £ Jr 

where dWt is a Brownian motion built on (S,M,P_b), Lt(x) is its local time and 
W(dx) is spatial white noise built on (Cl, JF, P). 

In the following, we follow their approach to show the convergence in distribution 
of £~ 1 ^ 2 (vi ie +V 2 , e )- To make the argument self-contained, we will provide all details 
and make appropriate modifications. 

To simplify the presentation, we will show the weak convergence of v\ yE /y/e and 
V 2 ,e/\f£ separately, and in the end it is easy to observe that the proofs combine to 
show the weak convergence of £~ 1/l2 (vi^ E + i> 2 , e ). 

4.1. A decomposition of the probability space. We decompose Cl as follows. 
Define 

(4.4) W e (x) := \fe<j)(x/E)/c, 

and since W e =$■ W in C(R), the family {W E } is tight. For any fixed <5 > 0, we can 
find a compact set K C C(R) such that for all s £ (0,1) 

P(W e G K) > 1 - S. 

Clearly K admits an open covering of the form 

u s eC(R){^ S C(R) : sup \h(x) - g(x)\ < 6}, 

xE[- S~ 1 ,S ~ 1 ] 

so we can extract finitely many deterministic function g \,..., gjy G C (R) such that 

K C U k=1 {h G C(R) : sup \h(x) — gk(x)\ < 6}. 

xG[- 6~ 1 ,S ~ 1 ] 

It is clear that we can further assume g & is bounded (the bound depends on S since 
K depends on S). Define 

(4.5) B k e = {uj G Cl : sup |V\4(x) — gk(x)\ < <S}, 

xG[- 5 ~ 1 ,<5 _1 ] 

and let B^' £ = B\’ £ and for any 2 < k < N, 

B s k ’ e = B S k e \ UjllBf 6 , 

so we have 

N 

P(uf =1 ^’ e )=^P(Bf’ E )>l-<5 

k =1 

for all e G (0,1). Let A s,s = Cl \ A k=1 B k e , so P(A <5,e ) < <5. Similarly, we define 
B k , B k and A s with W E replaced by W in (4.5). 
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The above decomposition of the probability space helps to “freeze” the random 
environment, i.e., with a high probability and a high precision, we can use finitely 
many deterministic functions to approximate W E . This helps to pass to the limit 
with only the “partial” expectation E#. 

4.2. Analysis of Vi >e . According to the decomposition of the probability space, we 
have 

Vi 


~r=- cE B {f'(x + M e t )y\> e (eX t/e 2 )} 
\je 


N 


= - KeAS.°cE B {f'(x + M^)W E (eX t/£ 2)} - ^ l ueB s, e cE B {f'(x + M t £ )W e (eA t/e 2)} 

k=l 

~(i) + (**)■ 

For (i), we first have 

\l^ A ^E B {f\x + Mt)W E {eX t/E .)}\ < W..E B {|W e (eX t/ea )|}. 

By taking E on both sides and applying Cauchy-Schwartz inequality and Lemma 3.2, 
we obtain 


(4.6) 


E{\l ueAS , E E B {f(x + Mt)W E {eX t/E .)}\} <^V{A^)EE B {\W E {eX t/e ^} 

<Vsti. 

For (ii), we write each summand as 

1 u6 b*-Eb{/ , (* + M^W E (eX t/E ,)} 

= ^ujeB 6 k ’ e ^ B {f'( X + Mt)yV E (£X t / E 2)(l\ E x t/e2 \>6- 1 + l|e.Y t/e 2 | <<5 1 )} • 

For the first part, we sum over k and write 

I E B {f'(x + Mt)W e (eX t / e 2 )l| £ x t/e2 j>( 5 -i}I ^ E B {|W e ( £ X t/e 2 )|l| £Xt/£2 |> 5 -i}, 
then the same proof as in Lemma 3.2 leads to 

(4.7) EE B {\W E (eX t/e *)\l leX < [ ^\\e~ c ^ 2/t dx. 

For the other part, we write 

1 cjeB 6 k ’ sE B{f'{x + M^)W E {£X t/E 2)l| e x t/e2 | < <5-i} 

=1 cjeB 6 k ’ eE B{f'(x + M t e )(W e - g k ){eX t / E 2 )l\ E x t/e2 \ <s -i} 

+ + XIt)g k (eX t / e 2) l| e jv t/E 21<5— 1 } 

:=(iii) + (iv). 

For (in), when u S B^’ £ , \W e (x) — gk(x)\ < 5 for |x| < <5 _1 , so 

( 4 -8) l w esf’ e |E B {f'(x + M^)(W e - g k )(eX t / e 2 )l\ e x t/e2 \<s- l }\ ^ <51^^. 

For (iv), we apply the quenched invariance principle of eX t / s 2 (and M t E ) to obtain 
that for almost every w, 

E B {f(x + M^)gk (eX t / s 2 ) 1 | e x t/s2 1 <5-i} 

-+E w {f'(x + dWt)gk(dWt) 1\ aWtlKS- 1 } 

as s —> 0. Furthermore, by the weak convergence of W 6 =>■ W, we have 


(4.9) 


1 


u<£Bl’ e 


u>eB\ 
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in distribution as e —> 0 , since the measure induced by W on C(R) does not sup¬ 
port on the boundary of the set {h £ C(R) : sup xe [_ 5 _i 5 -ii | h(x) — g k {x)\ < <5}. 
Therefore, we have 


1 ueB^ eE s{/'(a; + M^)g k (eX t/e 2 ) 1 | e _v t/e2 1 <( 5-i} 


in distribution as £ —> 0 . 

To summarize, we can write 


vi ,e{t,x) 

V~e 


N 


^ 1 ueB s k ’ ecE B{f'(x + Mt)gk{eX t / e 2)l\ eXt ^ 2 \ <5 -i} + R s>e 

k =1 


with the first part converges in distribution and EEs{|i? 5 , e |} — > 0 uniformly in e as 
S —> 0. Now if we write 

E w {f'{x + aW t )W(aW t )} =1 ueA sE w {f'(x + dW t )W(dW t )} 

N 

+ y i ueB^w {f'{x + aW t )W(aWt)l^ Wt \> s -i} 

k =1 
N 

+ y / 1 ojeB s k ^w{f(x + aW t )(W — g k )(&Wt)l\ 9Wt \ <s - 
fc =1 
N 

+ y^ 1 u£Bl^w{f(x + dW t )g(fc(d-Wt)l|CTVUt| < 5-i}, 

fc=1 

by the same discussion, we have similar estimates for the first three terms in the 
above expression as in (4.6), (4.7) and (4.8). Now we only need to send S —> 0 to 
conclude that 

(4.11) vi,eM ^ cE w {f( X + aW t )W{aW t )} 

\JE 

in distribution as e —> 0. 


4.3. Analysis of V 2 , £ - By the proof of Lemma 3.3, we have 
(M e ) t -at 3 f t/e2 


a2 


V~e 


=£ 2 


= - £ 2 


V(X s )ds 


ft/e 2 

${X t/e *)+el / ${X a )a{X a )dB„, 
Jo 


with 1 /> given by 


fe = -? r^ldy. 


a Jo a(y) 

Since W e (x) = y/e(j){x /e)/c, we can write 


(4.12) 


(M £ ) t -at _ 2c f [ eX */* 2 We{y)_ dy _ r 1 mjeXs/e 2 ) \ 


7,1 




%/e) io d(£X s/s 2 /£) 


where B s := £B s / e 2 . 

The idea is the same as for ui i£ , i.e., we decompose the probability f l as 

n = (u^^ e )u^ 6 

and for w £ B s k ' £ , we use the deterministic function g k to approximate W 6 and pass 
to the limit by the invariance principle of £X t / E 2 . 
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First, for w £ A S,E , by Cauchy-Schwartz inequality and Lemma 3.3, we have 

(4.13) E{|l weA ,,eE B {f"(x + MJ)((M E ) t - at)/V~e}\} < V® 

Secondly, we apply Cauchy-Schwartz inequality and (3.2) to derive that 

|E s {l suPae[o t] \ s x s/e2 \>s-if"(x + Ml){{M e ) t - at)/V~e}\ 2 
<E B {l suPse[o t] |EXs/£2| > 5 _ 1 }E B {((M e ) t - at) 2 /a) 

<e- c5_2 /*E B {m e )t ~ at) 2 /e}. 

By taking E on both sides and applying Lemma 3.3, we obtain 

(4.14) |EE B {l sup o \ eX \>s-if"(x + M E )((M E ) t - dt)/V~e}\ < 


Now we dehne 


G(W e ) := 


( M e ) t - at 


re-X, 


= 2 ca 


t/e2 W s (y) [* W e (sX s/e 2 ) 
;dy - 1 


dB s 


/0 d{y/e) J 0 cr(eX s/£ 2/e) 

and consider the error induced by replacing W e —> gk when w £ B ^ e . Since 

\We(y) ~ 9k{y)\ < 8 

for y £ [—(5 -1 ,<5 -1 ] and w £ B^’ e , we have 
(4.15) 

l ueB ;.«|EB{lsu P . 6[0 , t] I eX.,*\<6-'f'\x + M E )(G(W e ) - G(g k ))}\ 

^IE B {|eX t / e 2 |(5} + ^E B {j^ (W e - g k ) 2 {£X s / e 2)\\ eXa/ ^\ <s -ids} 

~ 1 wes*’ e {^BileXt/^S} + SVt^j <1 ^^eSVt. 

Here we have used the simple fact that 


, r* (W e -g k )(sX s/e2 ) 
i S u Pse[0 , t] | cX e/e2 \<6-A J o HeXs/e 2 /e) 


dB* 


< 


L 


tATs (W e - g k ){eX s/e2 ) 


dBs 


/ 0 d[EX s j e 2 j e) 

and Ito’s isometry, where := inf{s > 0 : eX s / e 2 > <5 -1 }. 

By combining (4.13), (4.14), and (4.15), it is clear that 

N 

, -sup sg[M k-Y s/e 2 |<<5 + Mt)G(g k )} + Rs,e 


(4.16) ^ = 1^,E fl {l. 


with Rg iC — )> 0 in L 1 )^.) as <5 —y 0 uniformly in e. 

Now we consider 

|eA' a/e2 |<<5- 1 / i x + )G(g k )} 

for fixed k. As before, for 1 => l^gs^ in distribution by the weak convergence 

of W £ =>■ W, so we only need to prove the convergence of 

E B {l SUPae[tM] |eA a/E 2 |< S-T-f'ix + M t)G{g k )}- 


G{g k ) =2ca 


r e X.,o 


t/e2 g k (y ) , [‘ gk{sX s/e 2 ) 


-dy - 


0 d(y/e) Jo a(£X s/E 2 /e) 
:=2cd(TZ E t - M E t ). 


dBs 


Let 
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For IZf , by ergodic theorem, we have 


1 


f eX t/e 2 


svPoejo.t] \zX a / s 2\<6 


9k{y) 


1 1' 

- - ) dy -> 0 


i(y/e) 


as £ —> 0 for almost every u and B s . Thus by the quenched convergence of sX t / e 2 
aW t in C(R+), we have 


(4.17) 


®s{lsup sg[0|t] \eX B/e 2 \<S~ 1 f ( x + Ms ffit} 

^ pd'Wt 

>E w {l sup \ 9Ws \ <s -if"{x + aW t )~ / gk(y)dy} 

a Jo 


for almost every w. 

For the continuous martingale 


Ml = 

its quadratic variation is given by 


/ 


9k i~X E / E 2 ) 


o d(sX s j e 2 /c) 

* 9k( eX s/eO 


dB. 


(M £ ) t = f 

Jo 


o d{sX s j e 2 /s') 


ds, 


and recall that in the proof of Proposition 3.1, sX t / e 2 is decomposed as sX t / e 2 = 
R E t + Mf with 


M E = a 


1 


dB s , 


1 0 d{sX s j e 2 / s) 
so the joint quadratic variation of M.\ and Mf is 

rt 9k(£X s / e 2) 


(M E ,M E ) t = a f 
Jo 


ds. 


1 0 d{sX s j e 2 / s) 

The following lemma shows the joint convergence of M E t and Mf. 
Lemma 4.1. For almost every ui, 


1 /'* 

[M e t ,M E ) => (- / g k (aW s )dW s ,aW t ) 
a Jo 


/ 0 

m C(R+) as e —> 0. 

Proof. We first consider the process on C(R + ) defined by 


h e {t) := [ 
Jo 


1 


0 d{sX s j e 2 J s) 


ds = s' 


ft/e 2 

"L 


0 a(T Xs u) 


ds. 


Since a^ 1 is bounded, h E {.) is tight in C(R+), and by ergodic theorem, its finite 
dimensional distribution converges to that of h(t) := t/a. Therefore, for almost 
every ui, we have h e => h in C(R+). 

Secondly, since sX t / e 2 =s* crW* in C(R + ), we have g\(sX s / e 2 ) =>• g k (oW s ) in 
C(R+). Now by [13, Lemma 3.5], the following mapping is continuous from C(R + ) x 
C(R+) -t C(R+): 

(hi, h 2 ) 1 —>• / hi(s)dh 2 (s), 

Jo 

when h 2 is increasing and C(R+) is equipped with the locally uniform topology. 
This implies 


(M E ) t = [ 
Jo 


9k( £ X s/ ^0 1 

1 ~ / v , . G?S _ 
0 X s / s' 2 / £) & 


g k {aW s )ds 
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C(R+). Similarly, 


(M £ ,M e ) t =d [ 9 f £ * s,e ^ ds =» [ g k (dW s )ds 

J 0 a \£X s / £ 2 /E) J g 


/q a(eX s / s 2 /e) J 0 

in C (M+). 

Now given the fact that 

(M e ) t => at, 

we conclude by martingale central limit theorem that 

(M e t ,M e t )^{- f g k (dW s )dW s ,aW t ) 

a Jo 

in C(M+) for almost every ui. The proof is complete. 

From Lemma 4.1 it is clear that 


□ 


1 


su Ps£[0,t] \z X s/e?\< 5 1 f ( X )M-t 


►1„ 


\<tw t \<8-if"(x + aW t )- f g k (aW s )dW ( 
a Jo 


su P S £[0,t] 

in distribution for almost every to. Now we only need to note 

r tAT6 9 k(£X s/e 2 ) 


Is 


-1-Mil < I / 

Jo 


~ ( y / \ dB s | , 

o 0 -(eX s / E 2 /e) 


Lsu Pae[0,t] \^ s/e i\<5 ~ 1 1 

which implies Es{l SU p 3e[o t] \ eX 2 |< 5 -i l-Mf | 2 } is uniformly bounded in e, so 

E B{lsup 3e[0 , t] \eX s/e2 \<8- 1 f"( x + M?)M e t } 

1 /■* 

tEw{l S up se[0it] \aW s \<6-if"( x + <rWt)-r J g k (aW s )dW s } 


(4.18) 


for almost every ui G Q. 

Combining (4.17) and (4.18), we obtain 
(4.19) 

N 

^weB^ e ® B {-*-sup 3eI o jt] h-Y s/e 2 1<5 -1 / (% + M t )G{g k )} 

fc =i 

/ raWt rt \ 

=^yi 1 weBf ]E tu{lsup 3e[0it] |stv S; |<5-i/"(a; + ^W ; t ) (2c g k {y)dy - 2cd g k (dW s )dW s \} 
in distribution as e —> 0. 

Now we send (5 —> 0 on the r.h.s. of (4.19) to obtain 

/ r&Wt rt \ 

r.h.s. of (4.19) — ► E w {f"(x + dW t ) ( 2c J W(y)dy — 2cd J W(aW s )dW s j } 

in L l (Gt). The discussion is the same as for (4.13), (4.14) and (4.15). 

To summarize, we have 


(4.20) 


V2 ,e(t,x) 


Kw{f" 


dWt 



<fW t 


W(y)dy — ca / W{crW s )dW s ]}. 
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Now we note that the proofs of (4.11) and (4.20) can actually be combined 
together to show that 
(4.21) 

Vl , e (t,x) +V 2 , e {t,x) 




- cE w {f'{x + aW t )W(aW t )} 


/ poWt rt \ 

+ E w {f"(x + aW t ) (c J W(y)dy — cd J W(dW s )dW s J }. 


By [13, Lemma 3.12], we have 


(4.22) 


r&Wt 


1 


W(y)dy — c<r / W(crW s )dW s = -R(0) 2 a 2 / L t (x)W(dx) 


with L t (x ) the local time of aWt . For simplicity we formally write the r.h.s. of the 
above expression as 

if?( 0)h 2 f L t {x)W{dx) = ^R(0)h 2 [ W(dW s )ds, 

2 J R Z J 0 

so (4.22) can be viewed as an application of Ito’s formula. 

Now (4.21) is rewritten as (recall that c = R(0)^a) 


(4.23) 


Vl ,e(t,x) + V2, e (t,x) 


- R(0)*aE w {f(x + dWt)W(dWt)} 


+ 


i 

-R(0)h 2 E w {f"(x + dW t ) / W(dW s )ds}. 
^ Jo 


5. Discussion 


5.1. The shift of the random environment and global weak convergence. 

The weak convergence in (4.23) is for fixed (t,x). Now we discuss the convergence 
of finite dimensional distributions of v(t : x) as a process in ( t,x ). 

From the beginning, we have shifted the random environment u> by t_ x / s1 and 
this is only used when applying martingale central limit theorem to obtain quenched 
weak convergence. The reason is that we need the environmental process u> s to be 
independent of e, and the shift of the environment enables the process co s = tx 3 x> 
to start from u instead of T x / e uj. 

Retracing the proof, the shift of the environment is used in proving Proposi¬ 
tion 3.1, (4.9), (4.17), Lemma 4.1 and (4.18) to get almost sure convergence, which 
is sufficient but unnecessary. For example, in (4.9), (4.17) and (4.18), a convergence 
in probability suffices to pass to the limit, which itself could come from an L 1 (fl) 
convergence, since the L 1 (f2) error does not depend on where the environmental 
process starts. Therefore, by almost the same proof, we have a convergence of fi¬ 
nite dimensional distributions, except that now the limit in (4.23) should encode 
the dependence on x. 

Without the shift, by the same proof (4.23) becomes 


(5.1) 

with 


Vl ,s(t,x) + V2,e(t,x) 

V~e 


=> v(t,x), 


v(t,x) = — R(0)^aE w {f(x + crWt)(yV(x + aW t ) — W(x))} 

I rt 

+-R(0)^a 2 E w {f"(x + dW t ) / W(x + dW s )ds}. 
^ Jo 


(5.2) 
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Comparing with (5.2) and (4.23), formally it is only to shift W(y) K > VV(a : + y), i.e., 

raWt 

' W(x + y)dy = W(x + aW t ) — W(x), 


raWt 

W(dW t ) = / W{y)dy 
Jo 


and 


W(aW s )ds H > / VV(x + aW s )ds. 


Since the spatial white noise W is stationary, the pointwise distribution is un¬ 
changed. 

Now we prove the “global” weak convergence, i.e., a spatial average of the ho¬ 
mogenization error u e — «hom of the form 

- 7 = [ (u s (t,x) - u hom (t,x))g{x)dx 

converges for every test function g £ C£°(M). 

First, we point out that the estimates derived in the proofs of Lemma 3.2 and 3.3 
grows polynomially with respect to |a:| when we have Xq = x/e instead of Xq = 0. 
Since g is fast-decaying, we have 

/ E{| (u e (t,x) - tthom(M) Vi t e(t,x) v 2 , e (t, x))g(x)\dx} <C yfl. 

Jr 

It still remains to analyze 

- 7 = [ (Vl,e(t,x) + v 2,e(t,x))g(x)dx. 

V £ J R 

The previous argument goes through except for some modification in the proof of 
(4.9), (4.17) and (4.18). Taking (4.9) for example, it suffices to show that 


(5.3) 


J g{x)E B {f'(x + Mf)g k (eX t/e 2 )l\ e x t/e2 \< 5 -i}dx 
-t / g(x)E w {f'(x + dW t )g k (x + dW t )l\ x+ ^ Wt \ <s -i}dx 


in probability as e —> 0. We emphasize that eX t / e 2 = x + + M t e and the 

environmental process w s = tx s oj starts from uiq = t x / £ uj. To prove (5.3), we 
consider the L 1 )!!) error 


(5.4) f | 5 (x)|E{|/i(a:) - I 2 (x)\}dx, 

Jr 

with 


h{x) = E B {f'{x + M^)g k (£X t/e 2 )l leXt/e2l<s -i}, 

h{x) = E w {f(x + dW t )g k (x + dWt)l\ x+ aWt\<5-i}- 


For every fixed x £ R, we can again shift the environment by r_ x / e without affecting 
the value of E{|/i(:r) — T 2 (a:) |}, and after the change of the environment, by the 
quenched invariance principle, E{|/i(a;) — T2 (m) |} — > 0 as £ — > 0. Since it is uniformly 
bounded, we only need to apply dominated convergence theorem to derive (5.4) —> 0, 
which further implies (5.3). 

To summarize, we have shown that 


/ \ U e (t,X ) Uhom(L*^) /. \ 

V e [t,X) = - 7 =- => v(t,x) 

V £ 

in the sense of Theorem 2.2, with v(t,x) given by (5.2). 
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5.2. PDE representations and a comparison with high dimensions. Now 

we discuss the individual terms coming from v(t, x). By (5.2), let us write 

v(t,x) = —R(0)^avi(t,x) + R(0)^av 2 (t,x) + -R(0)^a 2 v 3 (t,x), 

with 

Vi(t,x) = E w {f'(x + aW t )W(x + aW t )}, 

V 2 (t,x) = E w {f'(x + d-W t )}W{x), 

v 3 {t,x) = E w {f"(x + dWt) I W(x + aW s )ds}. 

Jo 

It is easy to see that v\[t,x) is the solution to the heat equation with a random 
initial condition, i.e., 

(5.5) d t vi = ^ad 2 vi, with iq(0,a;) = f(x)W(x), 
and 

(5.6) v 2 (t, x) = d x u hom (t, x)W(x). 

For v 3 , Lemma A.4 shows that it solves the SPDE with additive spatial white noise 
and zero initial condition: 

(5.7) d t v 3 = ^ad 2 v 3 + d 2 u hom (t,x)VV(x), with u 3 (0,x) = 0. 

We point out that v 2 corresponds to the first order fluctuation obtained by a 
formal two scale expansion. If we write 

u £ (t, x) = Uhom(t, x) + emit, x, -) + ..., 

£ 

then u\ (f, x, x/e) = d x Uhom{t, x)<f>(x/e) with (f> solving the corrector equation (3.3). 
Since y/e^x/e) scales to a two-sided Brownian motion (which is not stationary) 
when d = 1: 

y/e(j)( — ) => R(0)%aW(x), 
e 

we have 

£d x u hom (t,x)(j)(-) ~ y/£d x u hom (t,x)R(0)idW{x) = y/eR(0)^av 2 (t,x). 

£ 

The first order fluctuation given by v{t,x) is very different in high dimensions 
d > 3. Recall that u e — Whom ~ Ui, e + v 2 ^ e with v\ >e (t, x) = E s{/'(x + M t £ )i?f} and 

Rt = -e^X t/£ ,)+e4>(X „). 

When d > 3, we have a stationary zero-mean corrector [8, Corollary 1], so 

Eb{/'(x + Mt)£(f)(Xo)} ~ £d x u hoin (t,x)<f>(~), 

£ 

and 

|E B {/'(x + M £ )^(X t/e2 )}|« £ 

due to the fact that E{(/>} = 0 and the mixing induced by X t / e 2 when e is small. 
For v 2 , s (t,x) = ^E s{f"(x + Mt)((M e ) t — at)}, it turns out 

£~ 1 ((M E ) t - at) = a 2 £ / V(X s )ds 
Jo 

is an approximating martingale when d > 3, and is asymptotically independent of 
Mf. This implies \v 2tE (t,x)\ <C e. Combining these results, it was shown for fixed 
(f, x) that 

~ X 

U e {t , x) = Whom(t, x) + eVu hom (t, x) ■ </>(-) + o(e), 

£ 
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with o(e)/e —> 0 in L 1 (J7). Therefore, the pointwise first order fluctuation when 
d > 3 is given by eVrthom (t,x) ■ c/)(x/e), which only corresponds to V 2 (t,x) = 
dxUhom{t, x)W(x) when d = 1 . 

The following simple example illustrates the differences. Let u e (0,x) = £ • x for 
some fixed direction £ £ R d , so 

u e {t , x) = E s {£ • £X t/e i} = £ • x - eE s {£ • 4>(X t/s 2 )} + e£ • <t>(X 0 ) 

When a stationary corrector exists in d > 3, E#{£ • <j>(X t / e 2 )} —> 0 in L 2 (fl) as 
e —> 0 , and this is not the case by our proof when d = 1 . 

To summarize, the underlying diffusion process is so recurrent when d = 1 that 
the sample path is recorded in the asymptotic limit as e —> 0 , and all three terms 
in ui )£ + V 2 , s contribute to the first order fluctuation. When d > 3, we have 
sufficient mixing effects coming from the diffusion process, which leads to a different 
asymptotic behavior. 

5.3. An SPDE representation of v(t,x). At this point, our proof shows the 
limit v(t,x) is a superposition of three Gaussian processes V\,V 2 ,v 3 , and it turns 
out that they can be combined to form the solution to the SPDE given by (5.8): 

Proposition 5.1. v(t,x) solves 

(5.8) d t v = ^adlv - ^R{0)^a 2 d x (d x u iloin VV), with v( 0 , x) = 0 . 

We first give a heuristic derivation of (5.8). Recall that 

v 3 (t,x) = E w {f"(x + dW t ) I W(x + dW s )ds}, 

Jo 

and if we treat W as a function, an application of duality relation in Malliavin 
calculus shows that 

v 3 (t,x) = -E w {f'(x + uW t ) / W(x + aW s )dW s }. 

G Jo 

Furthermore, since v\(t,x) — V 2 (t,x) = E w{f'(x + dWt)(y\l(x + <f W t ) — W(x)), a 
formal application of Ito’s formula gives that 

1 R .. 

vi(t,x) - v 2 (t,x) - d 2 v 3 (t,x) = -d 2 E w {f'(x + dW t ) / W(x + dW s )ds), 

^ Jo 

so by recalling that a = d 2 , we obtain 

v(t, x) = — R(0)^a(vi(t, x) — V 2 (t, x ) — — av 3 (t , x)) 

= - ^(0)*a 2 (t > 3 + Vi), 

with 

Vi(t, x) = E\v{f(x + dW t ) f VV(x + dW s )ds}. 

Jo 

Since v 3 solves dtv 3 = ^dd 2 v 3 + d 2 u\ lom }V with zero initial data, the same argument 
should predict V 4 solves 

d t Vi = -ad 2 v 4 + 9 x Uh om VV, with u 4 ( 0 ,x) = 0 , 

hence v should satisfy 

d t v = ^ad 2 v - ^R( 0 )id 2 (d 2 u hom W + d x u horn W), with v( 0 ,x) = 0 , 
which leads to (5.8) if we write d x (d x U] mm \V) = S 2 WhomVV + S^MhomW- 
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The following is a rigorous proof of the above argument by introducing some 
mollification. 


Proof of Proposition 5.1. For fixed (t,x), the solution to (5.8) can be written as 

1 1 / pt 

•j(t,x) = - -R(0)i 


:= [ G(y)W(dy). 
JR 


q'a(t-s)( X - y)dy u h om (^5 y)ds J W(dy) 


It is straightforward to check that G £ L 2 (R) (since (t, x) is fixed, we have omitted 
the dependence of G on it). Define 


Ws{y) = [ h e (y - z)W(dz) 

Jr 


as a smooth mollification of VV. Here h e (x) = e 1 h(x/e) with h : R —> R + smooth, 
even, compactly supported and satisfying f R h(x)dx = 1. We can define 

v e (t,x) = [ G(y)W e (y)dy = [ G * h e (z)W{dz), 

Jr Jr 

and since G £ L 2 (R), G*h e —»■ G in L 2 (R), so v e —> v in L 2 (fl) as e — > 0. Now we 
show the L 2 (fl) limit of v e can also be written as a linear combination of vi,V 2 ,v 3 . 
First we rewrite v e as 


with 


S, x ) = - ^( 0 ) = 
:=(*) + (**)> 


(<) = -^m ^ 2 
(«) = -\m h d 2 


Qa(t-s)( x ~ y)d y (d y u hom {s,y)W £ (y))dyds 


<la{t-s)( x ~ y)dlu hom (s,y)W E (y)dyds, 


Qa(t- s )( x - y)d v u h om (&j y)W' e (y)dyds. 


It is clear that 


00 -t --R(Q)*a 2 v 3 (t,x) 


in L 2 (Q). For (ii), by the same proof as in Lemma A.4 we have 

I 1 pt 

(ii) = —- R(0)^a 2 E w {f'(x + aW t ) / W e {x + aW 8 )ds}, 
* Jo 


and an application of Ito’s formula gives 

rt r, r-x+SWt o P 

/ W' s (x + dW s )ds = - / W s (y)dy — - / W e (x + aW s )daW s . 

Jo a Jx a Jo 

For the second part, we apply the duality relation in Malliavin calculus and 
the fact that the Ito’s integral is a particular case of the Skorohod integral [18, 
Proposition 1.3.11] to obtain 

^w{f( x + dW t ) [ W e (x + dW s )ddW s } =aE w {f"(x + dW t ) f W e {x + aW s )ds} 

Jo Jo 

-dav 3 {t, x) 

in L 2 (D). 

For the first part, we write 

1 -x+aWt r 

Ws(y)dy= / (l[ x ,x+<rW t ]*h e )(y)W{dy), 
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so it is clear that 


E w {f'(x + aW t ) 




W E (y)dy} ->E w {f(x + aW t ) 
=vi(t,x) - v 2 {t,x) 


l[x,x+<rWt] {dy)} 


in L 2 (fl). The proof is complete. □ 

If we formally write in (5.8) that d x {d x u\ lom W) = <9 2 UhomVV + <9 x WhomVV, the 
term <9 2 MhomVV does not come from 1)3 since we have an opposite sign in (5.7). If we 
recall that v\ , v 2 comes from the remainder R% and V 3 comes from the martingale 
Mf, this indicates that the errors coming from the martingale decomposition need 
to be rearranged to obtain the correct representation given by (5.8). 


Appendix A. Technical lemmas 
Lemma A.l. E{|^(x)| 2 } < |x| and E{|'i/>(o;)| 2 } < |a;| 3 . 

Proof. Since 4>{x) = a J 0 ' V{y)dy and R(x) is the integrable covariance function of 
V, we have 

E {|0O C )| 2 } [ [ R(y - z)dydz < \x\. 

Jo Jo 

For ^(x), by (3.9) we have 

i>(x) = -- [ 4>{y)(V{y)+a~ 1 )dy = -2 [ (V(y) + a -1 ) [ V{z)dzdy 


W(z)| 2 } 


r x r x r y 1 r y 2 

< I i .-i 


l E {(^(yi) + a )(V(y 2 ) + a L ) v i z i) v ( z 2 )}\dz 1 dz 2 dy 1 dy 2 . 


/0 JO Jo Jo 


In the above expression, we need to control the second, third and fourth moments 
of V, which is a mean-zero stationary random field of finite range dependence. For 
the term with the second moment, we have 


fX fX pyi ry 2 


\R(zi - z 2 )\dz 1 dz 2 dy 1 dy 2 < \x\ 3 . 

Jo Jo Jo Jo 

The other cases are discussed in the same way by applying Lemma A.2. 
Lemma A.2 (Moment estimates). For any Xi £ R, i = 1, 2, 3,4, we have 


□ 


(A.l) 

and 


l E in^)}| < p{\xi ~ X 2 \) + p{\xi - x 3 |) + p{\x 2 - x 3 |) 


2=1 


(A.2) 


l E { II V(xi)}\ <p{\xi - x 2 \)p(\x 3 - xa\) + p(\xi - x 3 \)p(\x 2 - z 4 |) 
2=1 

+ p(|xi - Xi\)p(\x 2 - x 3 |) 

for some p : R+ — > R + satisfying p{r) < 1 A r~ p for any p > 0. 
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Proof. Since V is bounded, mean zero and of finite range dependence, (A.2) comes 
from [1, Lemma 3.1]. For (A.l), it is clear that there exists a compactly supported 
p : R + — > R + such that 

3 

mu V(xi )}| <p(min{|xi - x 2 |, |xi - x 3 |, |x 2 - x 3 |}) 

1=1 

<p{\xi - x 2 |) + p(\xi - x 3 |) + p{\x 2 - x 3 |). 

The proof is complete. □ 

Lemma A.3 (Estimates on local time). Let Lf(y) be the local time of a standard 
Brownian motion Wt starting from x up to t, then for any p > 1, 

/»00 

E{\LUy)\ P }<^ / qt{z)dz. 

J\y-x\ 


Proof First, Lf(y) has the same distribution as (y — x). By the strong Markov 
property of Brownian motion and distribution property of L® (0), we further have 

Lf (V ~ X ) ~ L®_T y _ x (0)l Ta _ x <J ~ M t -T y _ x ^Ty- X <t, 

where t v - x is the hitting time of another independent Brownian motion starting at 
zero and reaching at y — x, and M t is the maximum of Wt during [0,i]. Thus we 
have 

E{\L*(y)\P}= [ E{\M t - s \ p }p T v-*(s)ds<t% [ p T ^ (s)ds 
J 0 Jo 

= t%F(Ty- X < t), 

with p T y~ x the density of r y _ x . The reflection principle tells that 


y -x <t) = 2 


'\y-x\ 


qt(z)dz. 


The proof is complete. 


□ 


Lemma A.4 (SPDE representation). Letv(t,x) = E w{f{x+W t ) f 0 W(x+W s )ds}, 
then it solves 

(A.3) d t v(t, x) = -d x v(t, x) + u(t, x)W(x) 

with zero initial condition, and the function u solving dtu = \d x u with initial 
condition u( 0 ,x ) = /(x). 

Proof. The proof is similar to that of Proposition 5.1. First, we approximate the 
SPDE with a smooth equation. Then we use the probabilistic representation of the 
smooth equation and show its convergence. 

The solution to (A.3) can be written as 

v(t,x)= / / q t - s {x - y)u(s,y)W(dy)ds = / ( / q t - s {x - y)u{s, y)ds ) W{dy), 

Jo Jr Jr \J o / 

and we define v £ (t,x) as 

v e (t,x)= / / q t -s{x-y)u(s, y)W e (y)dyds = / ( / q t - s {x - y)u{s,y)ds\ W £ (y)dy, 

Jo Jr Jr \J o / 

with 

W e (y) = [ -h(- — -)W(dy) 

Js. e £ 

as a smooth mollification of W. It is clear that v e (t, x) -4 v(t, x) in L 2 (Ll) as e —> 0. 
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Since v e solves the equation 

d t V e = \ d l v e + 

by a probabilistic representation we can rewrite the solution as 

v e (t, x ) = Ew{ [ u(t — s,x + W s )W e (x + hF s )ds}. 

J o 

Since u solves the heat equation with initial condition u(0,x) = f(x), we obtain 

ve(t,x) =Eh/E s { f f(x + W s + B t _ s )W £ (x + W s )ds} 

Jo 

=E w {f(x + W t ) I W e (x + W s )ds} 

Jo 

=E w{f(x + W t ) [ W e (y)Lf(y)dy}, 

where L^{y) is the local time of x + Wt- 

By Lemma A.3, for any p > 1, E{|Lj(y)| p } can be bounded by some integrable 
function in y, and this helps to pass to the limit 

Ve(t,x) -> E w {f(x + W t ) f Lf(y)W{dy)} = E w {f(x + W t ) [ VV(a : + W s )ds} 

Jr Jo 

in L 2 (fi). The proof is complete. □ 
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